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We discuss real-space lattice models equivalent to gauge theories with a discrete non-Abelian 
gauge group. We construct the Hamiltonian formalism which is appropriate for their solid-state 
physics implementation and outline their basic properties. The unusual physics of these systems is 
due to local constraints on the degrees of freedom which are variables localized on the links of the 
lattice. We discuss two types of constraints that become equivalent after a duality transformation 
for Abelian theories but are qualitatively different for non-Abelian ones. We emphasize highly non- 
trivial topological properties of the excitations (fluxons and charges) in these non-Abelian discrete 
lattice gauge theories. We show that an implementation of these models may provide one with the 
realization of an ideal quantum computer, that is the computer that is protected from the noise 
and allows a full set of precise manipulations required for quantum computations. We suggest a few 
designs of the Josephson-junction arrays that provide the experimental implementations of these 
models and discuss the physical restrictions on the parameters of their junctions. 
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I. INTRODUCTION 

The ideal quantum computer 0, if it exists, is a 
very interesting object from a physical point of view. It 
is a system which may be prepared in any superposi- 
tion of 2 K (with K ~ 10 4 — 10 6 ) basic states during a 
very long time without decoherence. To give an order 
of magnitude, the decoherence should be much smaller 
than 10 -6 after a characteristic time scale of one opera- 
tion 0| . Furthermore, it should allow a set of controlled 
unitary operations with an accuracy of 10~ 6 or better. 
It is natural to think that in order to satisfy these strin- 
gent criteria one needs to have a physical system in which 
decoherence and errors generated by manipulations are 
exponentially small with respect to the system size or 
to any other control parameter. We shall refer to such 
properties as "ideal" in the following. The requirement of 
an ideal decoherence can be achieved in physical systems 
with topological order parameters 0, 0,11 which can be 
implemented in Josephson-junction arrays 0,0- How- 
ever, the physical arrays proposed so far do not allow the 
ideal implementation of the full set of operations Q, Q 
required for quantum computations. In fact, these arrays 
allow to perform ideally only the simplest operations: a 
flip of an individual bit and the change of the wave func- 
tion phase by ±1. It is crucial that they do not allow 
precise operations that mix the classical states, e.g., the 
Hadamard rotation Ujj that transforms the states |0), |1) 
into their mixtures (| 1) ± |0))/ y/2 requires an application 
of the correct external perturbation V for a given time t: 
Uh = cxp(— iVt) and thus is subjected to errors in both 
the values of t and V. 

Generally, an ideal operation involves an adiabatic 



change of some parameters that results in a generalized 
Berry's phase. Note that all alternative ways to operate 
the qubit are similar to a spin rotation in an external 
field. Such operations are not adiabatic and are con- 
trolled by the duration of the field pulse which can never 
be made precise. In fact, it is very difficult to imagine 
how one can possibly achieve the accuracy of the transfor- 
mation that is needed for quantum computation with this 
technique. In contrast, adiabatic transformations can be 
made, in principle, very precise. The main difficulty with 
them is to find a system in which the group formed by 
these transformation is big enough for a quantum compu- 
tation. The arrays discussed above are implementing the 
symmetry group H-i and thus it is not surprising that for 
them adiabatic changes of the external parameters allow 
only a very limited set of operations. One expects, fol- 
lowing the suggestion by Kitaev that one would get 
much larger sets of operations in systems allowing more 
complicated local symmetry transformations, in particu- 
lar, in the case of discrete non-Abelian groups [10(. 

The purpose of this paper is to consider lattice gauge 
theories built on more complicated groups than 7Li and 
to identify some that allow nontrivial adiabatic trans- 
formations. After having described the most promising 
ones, we discuss how they can be implemented, at least, 
in principle, in Josephson-junction arrays. All solid- 
state physics implementations of lattice gauge theories 
involve a lattice in real space but imply a continuous 
time. The natural framework to describe such models 
and their experimental realizations is the Hamiltonian 
formalism. This calls for an explicitly time-independent 
dynamics which can be achieved by choosing the scalar 
part of the gauge field Aq = 0. In this case, the vec- 
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tor part A^ describes the quantum state of the physical 
links. However, the Hamiltonian thus constructed still 
commutes with the large group of local time- independent 
gauge transformations. It is then convenient to use a de- 
scription which keeps track only of the underlying phys- 
ical degrees of freedom. This can be done by restricting 
the Hilbert space to the space of gauge-invariant states. 
Note that several studies have already pointed out the ex- 
istence of close connections between Josephson-junction 
arrays, lattice gauge theories and generalized statistics 
for elementary excitations. For instance, incompressible 
phases analogous to the fractional quantum Hall fluids 
(FQHF) of two-dimensional electronic systems have been 
predicted and studied in such arrays. One of them may 
be viewed as a FQHF of charges [Tlj induced by an ex- 
ternal magnetic field, whereas the other involves a FQHF 
of vortices |l2lll3| in the presence of a uniform fractional 
background charge on the superconducting islands. An 
appealing synthesis of these possible phases in the lan- 
guage of lattice gauge theories has been developed in Ref. 
|14j . As we shall show in detail, the lattice gauge models 
discussed here have the main advantage over the previous 
ones that effective interactions between elementary exci- 
tations are purely of statistical origin, and do not, for 
instance, depend on their precise locations. This is, of 
course, highly desirable in the context of quantum com- 
putation where the goal is to implement, with high accu- 
racy, a set of universal unitary transformations. 

This paper is organized as follows. In Sec. |H1 we 
review basic notions of continuous gauge groups and in- 
troduce basic notations. We emphasize the connection 
between the gauge-invariance constraint imposed on the 
physical states and the Gauss law constraint. Special fea- 
tures arising for non-Abelian groups are discussed. Sec- 
tion [ffi] focusses on the Abelian cyclic Z„ group. In this 
case, we show that imposing a local constraint on the 
product of group elements along the links of elementary 
plaquettes (kinematical Gauss law) is equivalent, after a 
duality transformation, to the requirement of local gauge 
invariance. We briefly discuss implementations of these 
two versions of Z„ models in Josephson-junction arrays 
and the existence of a topologically protected ground 
state when arrays exhibit macroscopic holes. Section llVl 
addresses the interesting case of theories associated with 
a finite non-Abelian group. In this situation, the usual 
local gauge-invariance constraint and the kinematical 
Gauss law constraint have to be treated separately, since 
no exact duality transformation is available to connect 
them. We show that both models present gapped low- 
energy excitations called fluxons which may be trapped 
in macroscopic holes of the arrays. We also discuss local 
violations of the constraints called pseudocharges, and 
show the existence of a non-Abelian Aharonov-Bohm ef- 
fect as a pseudocharge winds around a fluxon (and con- 
versely) . Section is dedicated to implementations of 
these two models, for the simplest non-Abelian groups 
such as the dihedral groups D3 and D±. These propos- 
als are based on the frustrated Josephson-junction arrays, 



where some basic elements have a 27r/n periodic potential 
energy as a function of the phase difference across them. 
The computations of the various energy scales driving 
the physical properties of these arrays, are presented in 
four appendixes. 



II. CONTINUOUS GAUGE THEORIES 

Let us start by a brief review of the properties of lat- 
tice gauge theories with continuous groups using the lan- 
guage that will be convenient for the discussion of dis- 
crete groups. Therefore, we consider a Hamiltonian H 
and a Hilbert space that are invariant under the local 
gauge transformations built from a continuous group G. 

More precisely, to any set of group elements {<?i}'s la- 
beled by lattice sites, we associate the transformation 
U {g} such that |*) -> U {g} \^) = and [H, U {g} ] = 0. 
For a given lattice, the gauge field A is attached to the 
links and takes its values in the group G. Local gauge 
transformations are then given by 
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In the familiar example of compact quantum electrody- 
namics (QED) A^ — exp(zajj), gi — exp(z^i) and gauge 
transformation in the variables (aij, 4>i) acquires a well- 
known form — » a^ + (pi — <j)j. Identifying quantum 
states I \P) with functions \& of classical field configura- 
tions {j4jj}, the condition of gauge invariance can be ex- 
plicitly written as 



(2) 



for any choice of g^s. In particular, we may choose trans- 
formations localized on a single site i, so that gj = e if 
j ^ i and gi = g, where e denotes the neutral element in 
the group G. Furthermore, since G is a continuous group, 
it is possible to consider infinitesimal g's. The constraint 
of local gauge invariance © can then be expressed as 
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where afj is the generator of the Lie group corresponding 
to the direction a in the vicinity of the neutral element 
of G. In the following, the notation stands for all 
sites j connected to i. The operator conjugated to afj is 
analogous to an electric field eg- = —id/dafj. In the basis 
where the electrical field e a is diagonal, Eq. © becomes 
a lattice equivalent of the usual Gauss law Ve Q = which 
also reads 
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where we have set Ef- = exp(iefA), t being an arbitrary 
real parameter. Notice that each site of the lattice corre- 
sponds to plaquette of the dual lattice, and the product 
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over the links emanating from a site becomes a product 
over the links encircling the corresponding plaquette of 
the dual lattice. This allows for a natural realization of 
conditions J2J) and Q in physical systems. 

At this stage, it is very important to point out a deep 
difference between Abelian and non-Abelian groups. In 
the former case, we may diagonalize simultaneously all 
the generators Efj since they mutually commute, which 
is impossible in the non-Abelian case. Because of this, we 
shall explore two different routes in this work. A natural 
possibility is to regard Eq. as the basic equation, 
and as we shall show, it is easy to use this construction 
for discrete non-Abelian groups. A second approach is 
to take Eq. 0} as the most fundamental. In this case, 
we assume that the physical states are also generated by 
classical configurations of group elements attached to 
the links, and we require these classical configurations to 
satisfy the kinematical constraint 

n^= e © 

Note that in this new scheme, and if G is non-Abelian, 
we can no longer interpret this modified Gauss law as a 
consequence of a local gauge symmetry. So, these two 
approaches yield essentially different theories, except in 
the Abelian case where their equivalence can be precisely 
established as explained in the following section. 

III. DISCRETE ABELIAN GAUGE THEORIES 



tween the various degenerate classical ground states dis- 
cussed above appear and usually lift completely this clas- 
sical degeneracy. Physical properties in the low-energy 
subspace spanned by these classical ground states are 
then well described by the following Hamiltonian: 

n-l 

H slt es = -Y,Y, r ( m ) T "( m ) = ( ? ) 
a£Ara=0 a 

where A denotes the physical lattice. In order to give a 
concrete illustration, we work on the geometry displayed 
in Fig. njwhere A is the hexagonal lattice and its dual A^ 
is the triangular one. Here, the operator T a (m) is defined 
by 

T a (m)\<p) a = \<p + 2Trm/n) a , (8) 

and r(m) is a tunneling amplitude which is expected to 
be maximal for m = ±l[n]. In this case, the Hamiltonian 





















/] \ 








.Uab=Uii 



Let us start with the usual Abelian case and con- 
sider more specifically a simple example where the dis- 
crete gauge group is the cyclic group 2„, which has 
also the major advantage that it can be easily realized 
in a physical system. Indeed, this can be implemented 
with Joscphson-junction elements for which the Joseph- 
son Hamiltonian is 2ir/n periodic in the phase difference 
{4> a — 4>b) where <p a denotes the phase of the supercon- 
ducting order parameter at site a. We discuss the ac- 
tual construction of such elements in Sec. [V] Due to 
this 27r/n-periodicity each link has n degenerate ground 
states characterized by an integer, u a b = 0, . . . , n — 1 
which is defined by 

A^ a6 = (j) a - (j) b = 2iru ab /n. (6) 

We would like to stress that this type of system exhibits 
a local ~K n symmetry defined by <f> a — > (f> a + 2Trm a /n with 
integers m a . Since we are dealing with quantum systems, 
it is also natural to interpret this symmetry in a basis 
which diagonalizes the local particle number operators 
conjugated to the phase operators. This symmetry then 
simply means that the particle number on each lattice 
site is conserved modulo n. An explicit realization for the 
n = 2 case has been recently proposed in Refs. 0, 0] • 

In real systems, charging energies induce quantum fluc- 
tuations of the phase variables. Tunneling processes be- 



FIG. 1: Schematics of the Josephson-junction array equiv- 
alent to the Abelian theory. Physical variables (supercon- 
ducting phases) are defined on the hexagonal lattice A: cj> a 
are phases of the islands while u a b describe the phase dif- 
ferences. Each bond represents a superconducting element 
which is 2n/n periodic in the phase difference <j> a — 4>b- The 
state of each bond is described by the discrete variable u a b 
according to Eq. ©. Alternatively, the state of the bond can 
be described by the variable Uij defined on the links of the 
dual (triangular) lattice Ad- 

H a is easily diagonalizable since its eigenstates are "plane 
waves" 

_^ n— 1 

\l) a = —f= / exp(i2irlm/n)\(j) + 2Trm/n) a , (9) 
Jn ^-^ 

v m=0 

and its eigenenergies are given by E(l) = — 2r cos(2irl / n) 
for I = 0, . . . , n — 1. Here, we have chosen tunneling am- 
plitudes r(m) equal to r for m = ±1 and otherwise. 
Thus the ground state \Q) = ® a \Q)a is given by a config- 
uration such that I = for all a. 

Although this model is very easy to solve since it only 
involves single-site operators which mutually commute, 
it is very fruitful to map it onto a genuine local gauge 
theory also based on the 7L n group. Indeed, this map- 
ping provides a convenient way to treat physical external 
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perturbations which are coupled to link variables. This 
construction is based on a description of physical states in 
terms of the link variables u ao . The crucial observation is 
that these variables are not independent as soon as there 
are loops on the lattice. More precisely, the sum of the 
integers w a b's along any oriented loop is constrained to 
vanish modulo n. A natural basis for our system is given 
by states of the form |\&) = ® a b\u) a b, where the u a b's 
satisfy the constraint. In the following, we shall refer to 
the description in terms of the <p a 's as "site" represen- 
tation, and to the latter (in terms of the u a b's) as the 
"link" representation. This condition on the u b's can be 
easily interpreted as the Gauss law for an Abelian gauge 
theory on the dual lattice. To show this, it is convenient 
to change the basis via the Fourier transform 



\p)ab 



1 
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e?q>(i2irpu/n)\u) ab , 



(10) 



M=0 



where p is an integer chosen modulo n. In this basis, a 
state \u) a b is an eigenvector for the translation operators 

^06 (?) w ith the eigenvalue exp(i27rp u/n) which provides 
a simple implementation of these operators in terms of 
the original fields. Thereafter, we will use the notation 
■V& (v) f° r the translation operator defined by 



T al ,(y)\x)ab = \x + y)ab, 



(11) 



for x = u,p. Note that [{T a ( f' '(p)] n = (np) = 1, so 

that the group composed of the s is also 2„. The 
constraint 

^ Uab = 0[n] <^ Y\_ cxp(i2iTUab/n) = 1, (12) 

plaq. plaq. 

for any elementary plaquette is equivalent to requiring 
the invariance of physical states under the action of the 
operator Yl plaq 7^(1) which shifts all the p a bS around 
the plaquette by 1. 

The last stage of the mapping onto a Z n gauge theory 
consists in expressing everything on the dual lattice. Any 
link ab on A separates two plaquettes which correspond 
to sites i and j of A^. The oriented links of A and A^ are 
thus in one to one correspondence. It is then possible to 
transform the p Q b's into the Pij's, and similarly the T^b's 
into the Ty's. We then recover the Gauss law written as 



(13) 



3(i) 



In other words, invariance of physical states under 

Ilpia<3 ^06 (1) i s equivalent to the requirement of invari- 
ance under the transformations pij — > pij + U — lj for 
any pairs of neighboring sites i and j belonging to A^, 
and any integers U and lj. This is exactly the equivalent 
of condition @ for a gauge theory associated with the 
group Z„ on the dual lattice where the gauge "fields" are 
the Pij's. 



The simplest Hamiltonian invariant under these trans- 
formations reads 



H = E E K«) II T ^ («) + E F ^ab 



(14) 



b(a) 



ab 



The first term gives quantum dynamics to the links and 
is equivalent to H S i te8 - The second term F(u ao ) de- 
scribes perturbations which lift the degeneracy between 
the classical ground states. In the physical implemen- 
tations in Josephson-junction arrays discussed in more 
detail below, it is usually due to the unprecise value 
of the flux through the elementary plaquette. For in- 
stance, in the realization 7] of 7Li theory it is equal to 
u a bEj5&/§o. Generally, this term breaks the invariance 
'Pa —>■ 4>a + 27rm a /n whereas it preserves the local gauge 
symmetry p^ — * j>y + /< — lj since it is diagonal in the 
| it) representation. 

The Hamiltonian H can also be written as 



H = -J2r(k 



n-l 



+ Pjk + Pk 



ij p=0 



(15) 

where f and F denote the Fourier transform defined by 



K(p) = exp(— i2irpu/n)K (u), 



(16) 



u=0 



for any function K. As usual, for a lattice gauge theory, 
the first term involves the "magnetic" flux (pij+Pjk+Pki) 
on the plaquettes ijk of A^, whereas the second term is 
off diagonal in the \p) representation and is associated 
with "electric" fields. 

In the absence of this latter term, the model can be 
solved and its ground state, in the \u) representation, 
reads 



\Q) = 



n s [ e ^ 

plaq. \(ab)£plaq. 



IlME 



n-l 



Ei 

«=0 



u=0 



It] i 



(18) 



where the first factor ensures the constraint l|12fl . The 
ground state can also be simply written in the \p) repre- 
sentation on the dual lattice 



where 



\g) = Y[Vi» ij \o)ij, 



p=0 j(i) 



(19) 



(20) 



is the projector on the gauge-invariant subspace defined 
by Eq. and where |0) stands for the states \p = 0). 
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The elementary excitations that we shall call fluxons 
correspond to plaquettes on Ad for which the correspond- 
ing flux is nonvanishing. Such an excitation is produced 
by an open string of operators 

/= [] ex P (i27ra o6 /n) = JJ T a ( fc p) (l), (21) 

(06)67 (ab)£7 

where the contour 7 goes over the links of the physical 
lattice A that is, over the centers of plaquettes of Ad on 
which the gauge field Pij is defined (see Fig. |2J|. As- 
suming r(m) = r if m = ±1 and otherwise, these ele- 
mentary excitations have energy E = 2r[l — cos(27r/n)]. 
Physically, in the system of Josephson junctions these ex- 
citations carry electrical charge 2e. We emphasize, how- 
ever, that in terms of the equivalent gauge theory these 
excitations carry no charge, they describe the fluctua- 
tions of the gauge field without matter. 




FIG. 2: The global topological invariant characterizing the 
ground state of the system is defined as a product of operators 
along contour which joins inner and outer boundaries (left 
bold line). Excitations (fluxons) correspond to the string of 
operators that are illustrated by the right bold line. 

Matter appears in this effective theory if one allows 
sites of the dual lattice where constraint is violated. 
In terms of the physical Josephson-junction system, it 
means that we have plaquettes of the physical lattice A 
where the sum of the superconducting phases is different 
from 0, that is vortices in the superconducting phase. 
Notice that due to 27r/n periodicity these vortices carry 
1/n flux quanta. Thus, it is natural to introduce sites i 
where ^Zj(i) = l[n]. In the gauge field language, it 
implies that gauge transformation at these sites changes 
the wave function of the gauge field according to 

II^GOI*) = exp(t27rp/n)|*>, (22) 



meaning that, at these sites, there is a pseudocharge 1/n. 
In order to create such an excitation above the ground 
state, one needs to apply operators 7y (1) along an open 
string on the dual lattice Ad. Because of the presence of 
this string, it is in general necessary to consider globally 
neutral charge configurations. This arises because the 
product of all the generators YljU) ^if over au sn - es * 
is the identity operator, or cquivalently because global 
gauge transformations do not modify the link variables 
Pij. When such a pseudocharge moves around a fluxon, 
it undergoes an Aharonov-Bohm effect with a dephasing 
exp(i27r/n). This is a natural result for a Josephson- 
junction array where a fluxon carries a physical charge 
2e and where a pseudocharge is a superconducting vortex 
with flux quanta ^o/n. The same result is obtained if 
we move a fluxon around a fixed pseudocharge as in the 
Aharonov-Casher effect. A more detailed discussion of 
this phenomenon is given in the context of non-Abelian 
gauge theories in Sec. IIVI 

The classification of the excited states as fluxons and 
pseudocharges is meaningful if these charges are rare. 
This is indeed the case if the energy needed to create 
a vortex of the superconducting order parameter (which 
becomes pseudocharge in the gauge theory language) is 
large compared with the Coulomb energy and also with 
the energy scale r of phase fluctuations. Consider qual- 
itatively the opposite case in which the tunneling am- 
plitude r is large compared to the potential energy of 
vortices and to their kinetic energy t. In a physical 
Josephson-junction array, this can be realized by con- 
structing the links of the lattice from rhombi [3 con- 
nected in series to an additional weak Josephson junction 
such that its Josephson energy ej and its charging energy 
tc satisfy ej < r and ej ^ ec 0- In these conditions 
all energy scales associated with weak junctions (ec, ej, 
vortex tunneling amplitude t) are small compared to r. 
Thus, in the leading approximation one can neglect the 
dynamics of the weak junctions and their contribution 
to the potential energy. The leading process remains the 
simultaneous flip of three adjacent rhombi described by 
the first term of the Hamiltonian i|14fl . The condition 
u ij = 0[ n ] no longer holds because weak links fa- 
vor the spontaneous excitation of vortices. Instead, it is 
replaced by the condition that all low-energy states have 
to minimize the kinetic energy part of the Hamiltonian 
(fH|) . i.e., to satisfy 

H^(m) = m. (23) 

b(a) 

This condition can be viewed as a Gauss law imposed on 
the sites of the physical lattice under the gauge group 
generated by T^f 1 [cf. Eq. ljT3|) ]. In these conditions, the 
physical array is expected to be in the insulating regime. 
It is therefore natural that the discrete version (|23[) of the 
QED gauge invariance © should be restored. Formally, 
the states satisfying the gauge-invariance condition i|23|) 
are described by the same theory as the states satisfying 



6 



Eq. I|13|) but the physical meaning of all quantities is 
reversed. Fluxons correspond now to the fractional vor- 
tices in the superconducting phase, while pseudocharges 
become the physical charges of the insulating array (mod- 
ulo n) . Finally, note that the degenerate states appearing 
in arrays with a nontrivial topology are distinguished by 
the product of 7y (1) operators taken over paths defined 
on the dual lattice. 

To conclude this section, we emphasize that when the 
gauge group is Abelian, there is a duality between the 
theory defined on the "real lattice" and the theory de- 
fined on its dual. This correspondence is established via 
the Fourier transform which relates the variables u and 
p defined on the real and on the dual lattices, respec- 
tively. In addition, when both structures are isomorphic, 
the gauge theory is self-dual, but it is not the case for an 
arbitrary system. In this sense, 7L n lattice gauge theories 
on an hexagonal lattice are not self-dual. However, the 
important characteristic required for the physical opera- 
tions discussed in our paper is not self-duality but duality 
which is always satisfied. 

IV. DISCRETE NON- ABELIAN GAUGE 
THEORIES 

In this section, we consider the case of a finite non- 
Abelian gauge group G. As already discussed in Sec. [HJ 
it then becomes necessary to consider models resulting 
from the gauge-invariance requirement and the kinemat- 
ical Gauss law constraint separately. 

A. Gauge-invariant model 

The requirement of local gauge invariance simply 
means that the wave function \^>) satisfies Us g \\ty) = \^>) 
for any configuration of the g^s, where U g is defined by 

UllA (24) 

Here, we have adopted multiplicative notation for the 
group law and we have introduced the left "translation" 
operator %j defined by 

T ij(h)\gh = \hg)ij- (25) 

As previously, link variables are oriented in the sense that 
\g)ij = \g The Hamiltonian invariant under those 

transformations has the form similar to Eq. I|15|) : 

H = - E r{9mk9ki) + E E P (9) T ii(g), (26) 

ijk ij g 

where r and F are now required to be invariant under 
inner group transformations g — > hgh^ 1 for any h in G. 
As before, f is associated with the energy cost of a fluxon- 
type excitation. Note that for a continuous group, the 



gij 's are naturally represented by matrices and f is often 
chosen as the trace function. By contrast, for a discrete 
group, we choose here the simplest case for which f(e) = 
r and otherwise, a choice that differs from previously 
where we had f(p) = 2r cos(2irp/n). As in the Abelian 
case, this model can be solved in the limit F — ► 0, where 
one can construct the ground state explicitly 

\Q) = l[ri®i j \e)i j , (27) 

i 

where 

is the projector on the gauge-invariant subspace. Note 
that this expression of the ground state is similar to the 
one given in Eq. (|19|l for the Abelian 7L n group. 

As in the Abelian case, the excitations are obtained by 
changing the value of the f on one plaquette. Neverthe- 
less, the non-Abelian nature of G prevents us to write 
down a compact form for the fluxon creation operator / 
which would generalize expression (|21ll . Instead, we first 
consider the state Ilffcz)^ Tki{g)®ij |e)y where the prod- 
uct goes over the links that are intersected by contour 7 
drawn on A that begins at the boundary of the lattice 
and ends at the excited plaquette. Then, one projects it 
onto the allowed Hilbert space getting 

l/) = n^ II ?«(s)®tf|e>«. ( 29 ) 

i {kl)e~i 

Since the projection operator JX;^ is gauge-invariant, 
one can insert any gauge transformation operator be- 
tween it and IT(jM)ey Tkl(jj) without changing |/). In- 
serting the global gauge transformation with a group el- 
ement h, we transform the string of Tki(g) operators into 
the string of Tki(g') with g' = h~ 1 gh. Thus, all group 
elements that belong to the same conjugacy class result 
in the same excitation with energy E = r. 

It is useful to give a more intrinsic description of these 
fluxons by considering generalized magnetic fluxes along 
closed paths. For such a path 7 C , the generalized mag- 
netic flux is defined as $(7c) — Yluj)etc 9*3' ^ n ^ e 
sence of fluxon, one has <fr(7 c ) = e independently of j c . 
In the fluxon state |/), one has <&(7 C ) = e only if 7 C does 
not encircle the plaquette on which the excitation is lo- 
calized. If the path encloses this plaquette once, then 
3>(7c) belongs to the conjugacy class of g, where g is the 
same group element as in Eq. I|29|l provided the orienta- 
tion of 7 C coincides with the one of the string operator 
used to build |/). Let us mention three important prop- 
erties of these generalized fluxes. For this purpose, we 
need to consider closed paths j c which begin and end at 
a fixed arbitrary site O. In this case, we shall say that O 
is the origin of j c or that j c is based at site O. We have 
then the following 
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1. If 7c and j' c are both based at site O, and if one 
may deform j c into j' c without going through any 
localized fluxon, then &(-f c ) — ®(l' c )- 

2. If 7 C and j' c involve the same closed loop but with 
different origins O and O' , then $(7 C ) and $(7 C ) 
belong to the same conjugacy class. 

3. < i ) (7c) is not gauge-invariant in general, but its con- 
jugacy class is gauge- invariant. 

For two or more fluxons, the quantum state of the pure 
gauge theory is no longer determined by the knowledge 
of the conjugacy class associated with each individual 
fluxon. This arises because if g\ is conjugated to g[ and 
92 to g' 2 , then gig 2 is not always conjugated to g'ig' 2 . Al- 
though all these conclusions are valid for an arbitrary 
group G, we shall explicitly consider here only the sim- 
plest non-Abelian group D n (n > 3) which contains n 
rotations ( k (fc = 0, . . . , n — 1) and n rotations combined 
with reflections tC, ■ The two generators ( and r satisfy 
the three relations: = e, t 2 = e, and t( = (~ 1 t, where 
e is the unit element of the group. For simplicity, let us 
choose the group D$ which contains six elements and 
three conjugacy classes: {e}, {C,C 2 }> an d {t, t(, t£ 2 }. If 
we take gi = g[ = gi = t and g' 2 = r£, then g±g2 = e and 
9W2 — C belong to different conjugacy classes. A com- 
plete description of the two-fluxon state requires then to 
specify the conjugacy classes of the fluxes associated with 
three sets of paths on A^: the paths which encircle each 
fluxon only once without going around the other, and the 
paths which encircle both fluxons only once (see Fig. |3J). 




FIG. 3: The closed paths 71 and 72 encircle only one fluxon, 
whereas 7 encircles both. The conjugacy classes of $(71), 
$(72), and $(7) are gauge-invariant and are also invariant 
under "smooth" deformation of the paths. However, $(7) is 
not conjugated to c I > (7i)3 > (72), except if all the paths have the 
same origin. 

Here, we have used the properties Q and @ of gen- 
eral fluxes given above. The non-Abelian nature of G 



induces then a notion of nonlocality when we combine 
several fluxons. In this sense, these objects do not behave 
as ordinary independent quasiparticles when we consider 
their global quantum statistical properties [Tof. 




FIG. 4: The path 72 is equivalent to 7^7271 which is not in 
the same class as 72 as this figure illustrates. Consequently, 
the fluxes $(72) and $(72) belong to the same conjugacy 
class, but they are in general different, in spite of the fact 
that 72 and 72 have the same origin O, that they both wind 
exactly once counterclockwise around fluxon 2, and do not 
wind around fluxon 1. 

More precisely, the state of the system is character- 
ized by the values of flux calculated for a set of closed 
paths that begin and end at a given site O of the lat- 
tice and avoid the fluxon positions. The value of flux 
is the same for all topologically equivalent paths (with 
the same origin), i.e., for the paths that can be smoothly 
deformed into each other without going over the fluxon 
positions. The set of topologically nonequivalent paths 
has a natural group structure (provided by path con- 
catenation), called the fundamental group of the lattice 
with n plaquettes removed, we shall denote it by 7Ti(n). 
As is well known |16|. 7Ti(n) is generated by the equiv- 
alence classes of n paths, each of which winds exactly 
once around only one particular fluxon. But for n > 2, 
this group is non-Abelian as illustrated in Fig. 0] No- 
tice that specifying a set of winding numbers for a path 
with respect to each fluxon does not determine the to- 
tal flux along this path which depends on the order in 
which the path encircles the fluxons. Formally, any con- 
figuration of link variables gij with n fluxons establishes 
a (non-Abelian) group homomorphism from 7Ti(n) to G, 
which sends the equivalence class of a path 7 based at 
O into $(7). Clearly, $(7172) = $(71) $(72). In the 
simpler case of Abelian G, this subtle effect due to the 
non-Abelian nature of tt\ (n) is washed out and fluxes do 
depend only on the set of winding numbers and not on 
their actual sequence along a path. Finally, we wish to 
characterize gauge-invariant states with n fluxons. As 
already noted, a gauge transformation g'^ = higijhj 1 



8 



transforms fluxes $(7) into $'(7) = ho^^hy 1 . Thus, 
equivalence classes of group homomorphisms from n\ (n) 
to G are in one to one correspondence with the gauge- 
invariant states for n fluxons. 

Let us illustrate this for n = 2 and G = D3. If flux- 
ons were described by group elements, we would expect 
6 2 = 36 states. But this clearly does not take the gauge- 
invariance condition into account. If we naively associate 
conjugacy classes to fluxons, this yields 3 2 = 9 states for 
two fluxons. It turns out that the correct counting is 11. 
We list below the corresponding equivalence classes of 
group homomorphisms from 7Ti(2) to D3. Since 7Ti(2) has 
two independent generators 71 and 72 (see Fig. 0J, this 
homomorphism is completely determined by the knowl- 
edge of $(71) and $(72)- In these notations, the eleven 
possible states of two fluxons are: 



$00 


= {(e,e)}, 


$01 


= {(e,C),(e,C 2 )}, 
= {(C,e),(C 2 ,e)}, 


$10 


$02 


= {( e ,T),(e,TC),(e,rC 2 )}, 


$20 


= {(T,e),(TC,e),(rC 2 ,e)}, 

= {(C,C),(C 2 ,C 2 )}, 
= {(C,C 2 ),(C 2 ,C)}, 


$11 


$11 


$22 


= {(t,t),(tC,tC),(tC 2 ,tC 2 )}, 

= {(r, rC), (r, r^), «, r), (r£ rC 2 ), (K 2 , r), 


$22 


(rC 2 





$ 12 = {(C, r), (C, rC), (C, rC 2 ), (C 2 , r), (C 2 , <), (C 2 , K 2 )}, 
$21 = {(r, C), K, C), (K 2 , C), (r, C 2 ), «, C 2 ), « 2 , C 2 )}- 

(30) 

Note that in this list there are two pairs of different 
two-fluxon states, $n, and $22, $22 which are 

produced by the combination of nondistinguishable one- 
fluxon states, namely, states with fluxes $1 = {£, £ 2 } and 
cf> 2 = {r,TC,rC 2 }. 

As usual, including static local charges in such a model 
requires to consider sites which transform differently un- 
der the gauge group. More precisely, these sites can be 
classified according to the representations of the gauge 
group. For a non-Abelian group, most irreducible repre- 
sentations are no longer one dimensional. This leads us to 
attach a finite-dimensional representation space spanned 
by basis vectors \a) to each of these sites. The corre- 
sponding representation of G in this subspace is described 
as usual by a set of matrices D a p(g). Now, let us assume 
that we have only one such site. The global wave function 
may be written as \^f) = J2 a \ a ) ® In this slightly 

generalized setting, the local gauge-invariance constraint 
becomes \^> a ) = J2p D a p (g)U g ), if the group element 
g is attached to a site carrying a nontrivial representa- 
tion. For the D n group that we consider here, the most 
interesting cases are obtained from the two-dimensional 
irreducible representation of the group. The wave func- 
tion of the ground state satisfying this constraint at one 
site o, occupied by a pseudocharge in the D(g) represen- 



tation may be written explicitly as 

\Q) = (il^j V \a) ®y |e)y, (31) 

where \a) is any vector in the charge Hilbert space, and 

As for the Abelian case, it is interesting to discuss the 
situation where a pseudocharge moves around a fluxon 
and conversely. The state combining a pseudocharge \a) 
located at site i and the gauge degrees of freedom in 
the state |\&) will be denoted by \a)i ® In a trivial 
background where \$>) = ®ij\e)ij, it is natural to assume 
that the charge hopping Hamiltonian H t acts according 
to 

ff t (|a)i® |*>) = -t £>>i®l*>. (32) 

m 

where t is a positive hopping amplitude. Requiring in- 
variance of H t under local gauge transformations com- 
pletely determines the action of H t for a nontrivial gauge 
field background 

H t (\a)i ® I*)) = -t'^2[D(g^ j 1 )\a)]j (g) |*). (33) 

m 

Because of the string involved to build the fluxon |/), 
it is clear that a pseudocharge moving around it expe- 
riences a generalized Aharonov-Bohm effect which is, in 
general, nondiagonal in the Hilbert space associated with 
the internal degrees of freedom of the charge. Conversely, 
we may consider the case of a fluxon associated with the 
conjugacy class of g and moving around a fixed charge. 
When the fluxon moves from one plaquette to an adja- 
cent one, we obtain the final state by modifying only the 
field variable gij on the link (ij) that has been crossed by 
the fluxon. This defines an operator Mij which is equal 
to Tij(g') where g' is a group element belonging to the 
same conjugacy class as g. Note that the choice of this 
representative element inside the class of g depends on 
the actual field configuration. Consider, for instance, the 
example of two plaquettes shown in Fig. [5] 



4 4 




1 1 



FIG. 5: Fluxon moving from plaquette (123) to plaquette 
(234) involves a modification of t/23- In the initial state (left), 
923^' = 924343, and in the final state (right) g^ n ' = (721313 • 
The operator which transforms the former into the latter state 

is M23- 
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Thanks to this construction, Mij commutes with all 
the local gauge transformations operators C/{/j} acting on 
the gauge degrees of freedom. Let us take an initial state 

V(\a)i ® |*)), where V = (U^ V?j V a . We move the 

fluxon along a closed path 7 C which encircles the origin 
o where the pseudocharge is located. This leads to the 
state: 



n m h Ni«>«>®i*» = ( 34 ) 

V{\a) ® I J] MiA |*». 

Since \^) and M.ij\^>) both describe the same fluxon, 
there is a gauge transformation U^y connecting these 
two states. We choose {h} such that hi — e for i outside 
the area enclosed by the fluxon trajectory. In this case, 
h a is conjugated to the flux $ = g associated with the 
fluxon. Our final state is then V[D(h- l )\a) ® |*)]. We 
thus recover the same non-Abelian transformation after 
a closed cycle as for the moving charge situation. 

Furthermore, non-Abelian Aharonov-Bohm effects also 
appear when a fluxon goes around another fluxon |l7j . 
This situation is shown in Fig. where fluxon 1 car- 
rying initially the flux $1=51 goes around fluxon 2 
(with initial flux $2 = .92) once counterclockwise. These 
fluxes are defined as a product of group variable around 
contours 71 and 72 starting and ending at point O or 
any other contours topologically equivalent to them. The 
state with these fluxes can be represented by the configu- 
ration of gij that contains two strings against the "flat" , 
gij = e, background. As the first fluxon moves around 
the second one, its string eventually crosses the string of 
the static fluxon. After the crossing the string of the first 
fluxon changes into the string of G\ = <7 2 ^ 1 3i52 elements. 
The new fluxes $1 and $2 become $1 = g 2 X 9\92 and 
$2 = 92 1 9i X 929W2- Since g^g^n = 9 2 X 9\ X 9\9\9i, this 
state can be reduced to the initial state by a global gauge 
transformation with the group element g = gig 2 , if there 
are only two fluxons in the system. But if some other 
fluxons are present somewhere else in the system, besides 
the two already discussed, they might transform nontriv- 
ially under the gauge transformation with g: <E> ^ g _1 <^g 
and in a general case the final state cannot be reduced to 
the initial state by the gauge transformation. This pro- 
vides another example of nondiagonal Aharonov-Bohm 
effect, even in the absence of external charges. These 
considerations have been generalized to composite ob- 
jects of fluxons and pseudocharges usually called dyons 
|10| . The appropriate mathematical framework has been 
developped in terms of representation theory of the so- 
called quantum double algebra associated with the group 
Il9l . We shall not elaborate further on these inter- 
esting topics here. 




FIG. 6: The effect of motion of one fluxon around another. 
Two fluxons, 1 and 2 carry initial fluxes $1 = g\ and $2 = 
<?2 that are calculated as a product of group elements along 
contours 71 and 72 beginning and ending at the same point O 
of the lattice. These fluxons can be represented as strings of g\ 
and 32 (full thick bonds) against the background of unit group 
elements. When fluxon 1 moves around fluxon 2, it leaves 
behind a string of group elements (dashed and dotted bonds) : 
g\ initially, until it crosses the g2 string, and afterwards G\ = 
92 1 9 1 9 2 - After this circular motion the fluxes become $1 = 
g 2 ~ 1 gig 2 and $2 = 32~ 1 fi~ 1 323i5'2 respectively. 



As for Abelian models, the systems with one or sev- 
eral holes have a degenerate ground state. In the limit 
of zero off-diagonal term in Hamiltonian (|26[) the fluxons 
are trapped forever in the holes and have zero energy. 
The degeneracy of the ground state with one hole is thus 
equal to the number of different fluxons, i.e., to the num- 
ber of conjugacy classes. The operator that distinguishes 
these ground states is Yic 9 wnere contour C encirles 
the hole. A small but nonzero off-diagonal term in the 
Hamiltonian leads to an exponentially small amplitude 
for the fluxon escape from the hole. The operator lead- 
ing to this escape process is the product of T a b(g) along 
the contour 7 on lattice A that connects inner and outer 
boundaries as shown in Fig. [21 This is the same opera- 
tor that moves fluxons around the system, as discussed 
above. The gap for the fluxons implies that this process 
is only virtual and its amplitude is exponentially small. 
Therefore, the resulting degeneracy splitting is also expo- 
nentially small. Note that (as discussed above) the choice 
of the group element g in fluxon motion operator 7~ a t,(g) 
depends on the link variables. The contour 7 crosses 
once any contour C and thus changes the flux trapped 
inside the hole. Since these fluxons undergo non-Abelian 
Aharonov-Bohm effects, they may be used to construct 
a complete set of elementary gates for quantum comput- 
ing, provided the group G is "sufficiently" non-Abelian. 
This has been recently analyzed in great detail by Mo- 
chon, first restricting the discussion to states involving 
only fluxons [2(j, and then by extending it to the case 
where pseudocharges are also involved |21|. This latter 
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construction is specially interesting for our purpose, since 
it shows that the small D 3 group (whose practical imple- 
mentations are discussed in Sec. EJ is already capable of 
universal quantum computation. 



B. Kinematical Gauss law model 

As explained at the end of Sec. [HJ a different model 
may be constructed by assuming a modified Gauss law 
in which the constraint is diagonal with respect to the 
variables a,-,-, ft reads 



(35) 



and can be directly compared to Eq. |JSj| . This constraint 
allows for a simple interpretation which is motivated by 
the physical systems discussed in Sees. IIIII and IVl First, 
we notice that a constraint associated with the sites of 
Arf is also naturally attached to the plaquettes on A. We 
may therefore also write it as Ilp/ag 9 ab ~ e - Second, the 
constraint that the product of group elements around the 
plaquette is equal to e implies that one can define group 
elements on the sites such that g a b — 9a% X ■ Further, 
these group elements are defined up to a global trans- 
formation g a — ► g a h. For the Abelian systems discussed 
in Sec. IIIII the group elements on the sites of A had a 
transparent physical meaning: they corresponded to the 
absolute values of the phase of the superconducting is- 
lands. 

The simplest Hamiltonian preserving the constraint 
is 



H = £ £ <9) II T M + £ F (9ab 



(36) 



6(a) 



ab 



Note the close similarity with the Abelian model defined 
in (I14|l . However, for a general non- Abelian group, we 
do not know how to reexpress the Hamiltonian i|36|) in 
terms of a local gauge theory on Ad- Indeed, the key tool 
in this mapping for the 2„ group was the existence of 
the Fourier transform for functions defined on 2„. Here, 
we do not see any obvious way to generalize this trick 
to functions defined on a non- Abelian group. Neverthe- 
less, if F = 0, the Hamiltonian (|36|l does exhibit a local 
symmetry based on the sites of A: g a b — > h a g a bh~^ x pro- 
vided r{g) depends only on the conjugacy class of g. In 
this case, the model can be solved for F — > because, in 
this limit, we may interpret the Hamiltonian as a sum of 
mutually commuting terms attached to the sites of A 

H sltes =-Y,Y, r (9)Ta(g)=Y, H - ( 37 ) 



Solution for each site depends somewhat on the structure 
of the coefficients r(g). For the simplest case of positive 
r(g), the ground state is simply \G) a = J2 g lff)a- On the 
whole lattice, the ground state is thus given by: \Q) = 



<&a\G)a- This ground state can be also represented in 
terms of the original link variables: 
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iim n 

plaq. \(ab)£plaq. 



g a b - e 



9ij 



£i»>* 



(39) 



We underline that these expressions are the generaliza- 
tions of l|17fl and (|18fl for an arbitrary non- Abelian group. 

This model exhibits local excitations attached to the 
sites of A or equivalently to the plaquettes of A^. They 
are analogous to the Abelian fluxons constructed in 
Eq. i|21|) where the site Hamiltonian H a in a system 
of Josephson junctions involves charging energies. The 
Hilbert space attached to the site a forms the so-called 
regular representation of G where the group element h 
is represented by the left multiplication operator T a (h). 
As already pointed out, these operators commute with 
H a so that the eigenstates of H a belong to stable sub- 
spaces for the regular representation of G. For exam- 
ple, if we choose G = D3, its regular representation is 
six dimensional and decomposes into two different one- 
dimensional representations and two isomorphic copies of 
the two-dimensional representation of D3. At this stage, 
we should notice that H a also commutes with the rep- 
resentation of G by right multiplications (the element h 
being represented by lZ(h) so that lZ(h)\g) — \gli~ 1 )). 
Because of this extra symmetry, the spectrum of H a is 
composed of two singlets and one quadruplet. 

In fact, it is possible to generate the whole set of eigen- 
states of H when F = 0. Let us choose a reference state 
\Q) in which the constraint Ilpiaq 9 ab = e i s satisfied for 
all plaquettes. Now, for each site, we consider an eigen- 
function 4> a (g) of H a with energy E a . More precisely, we 
require 



E a ci> a {g) = -Y,r{h)(t )a {h- 1 g). 



(40) 



We can then construct an eigenstate of H with the 
total energy E = J2 a E a as follows: 



i*)=riE^(ff)ii r ^(-9)i fi )- 

a g 6(a) 



(41) 



In the case where (j)a(g) is independent of <?, we recover 
the ground state \Q) given by Eq. (J^HJ- Introducing the 
gauge transformation operators U{ g y as before (where the 
<7's now live on the sites of A), we may also write 



\*) = HY,M9a)u {g} \n). 

a {9} 



(42) 



This shows that using a gauge transformed reference 
state instead of |f2) is equivalent to replace the 



11 



site wave functions 4> a {g) by 4> a (gh~ 1 ). As already dis- 
cussed, acting on the tf> a 's by these right multiplications 
produces eigenstates with the same energy. Only sites a 
where <j) a belongs to a multidimensional representation 
of G generates new states when a gauge transformation 
with h a 7^ e acts on 

This language is specially useful to describe a differ- 
ent type of excitation, analogous to the pseudocharge 
in the gauge-invariant model. They are produced in a 
physical system in which the constraint is modified on a 
given plaquette so that JlpZag 9ab — u 011 this plaquette. 
Note that as it stands, this constraint is not preserved 
by the action of the Hamiltonian (|36|l . but the conjugacy 
class of Y\ plaq g a b is preserved for any plaquette. This 
modified constraint may be incorporated in the general 
construction presented above provided we change the ref- 
erence state. For instance, we may consider a string 7 on 
h.d which runs from the center of the excited plaquette 
to the boundary of the system. Our previous state |f2) is 
then replaced by ri(ab)e7 T a b(u)\Q) ■ Note that in the con- 
text of physical implementations, as already mentioned 
in Sec. IIIII it is more appropriate to call these excitations 
vortices. 

As in the gauge-invariant model, it is interesting to 
study the appearance of nondiagonal adiabatic transport 
of a degenerate fluxon excitation around a pseudocharge. 
As before, let us describe the situation according to both 
viewpoints where either one excitation is fixed or the 
other. For a mobile pseudocharge, we have to transform 
adiabatically the reference state. Starting from 

l*)m=nE^W^} II r ab(u)\n)), (43) 
a {g} \(afc)67 / 

we obtain after one closed path f3 of the pseudocharge 
around the fluxon: 

\*)fin=HJ2M9a)U {g} I J] r ab(u)\n)\. 

a {g} \(ab)G 7 +/3 / 

(44) 

The projection operator fl a Y,{ g } 4>a{ga)U{ g } ensures 
that the resulting state has exactly the same energy as 
the initial one, so that the motion of pseudocharge does 
not excite fluxons. To make a closer analogy with the 
gauge-invariant model of Sec. II V Al we note that this 
wave function can be also written as 

\(ab)ef) J 

where operators Arab's that move the pseudocharge are 
defined on the links of A exactly as the A^y's were de- 
fined on the links of A^ in Sec. IIV Al these operators 
have the same effect as T ao {u) on the flat (|f2)) state and 
they commute with Hb( a ) ^ab(g) for any site a and any 
group element g. Using the same reasonings as for the 



gauge-invariant model, we see that ri(ab)e,3 M a b acts in 
the same way as a gauge transformation U^} on the state 
Since h has to be nontrivial (more precisely conju- 
gated to u) at the site where the fluxon is located, \ty) fi n 
is degenerate with \^)i n , but these two states may be 
linearly independent if the fluxon is created in a multidi- 
mensional representation of G. 

Alternatively, consider a motion of the fluxon which is 
due to a small hopping term, t(g a b) [arising, e.g. from F 
in Eq. JSSJ]- Let us denote by (j>f(g), the single-site wave 
function associated with the fluxon. In the absence of a 
pseudocharge, the fluxon state at site a corresponds to 

\*)a=Y,M9a)U {g] \n). (46) 

{g} 

Acting on this state with the operator t(g a b) produces 
the state |*) 6 provided ^(fl-a)^^ 1 ) = 4>f{gb)- H2- 
In the presence of a pseudocharge, is replaced by 
Il(ab)G7 1~ab(u)\£l) . If the link (ab) is on the string 7, one 
has: 

{9} \(ab)e 7 / 

(47) 

So, as the fluxon hops from a to &, its internal wave func- 
tion is modified by a right multiplication associated with 
an element conjugated to u and there is a complete agree- 
ment between both descriptions. 

Note the close similarity between these results and the 
ones for the gauge theory model, provided we exchange 
the roles of pseudocharges and fluxons. For instance, a 
pseudocharge in the kinematical Gauss law model expe- 
riences a nondiagonal Aharonov-Bohm effect when an- 
other pseudocharge is moved around it, even in the ab- 
sence of fluxons. Notice that although the pseudocharges 
in the Gauge theory are usually associated with inde- 
pendent degrees of freedom, this distinction (with flux- 
ons in kinematical model) is mostly formal. It com- 
pletely disappears in the alternative (albeit less popu- 
lar) gauge theory formalism in which one keeps only the 
dynamical variables gij on the links and creates pseu- 
docharges by replacing the projector V Q in Eq. i|31|) by 
Eh X{h) Ylj(o) 1~oj{h) where \{h) is the character of the 
representation of G associated with the pseudocharge at 
charge o. This approach emphasizes the analogy between 
the two situations and is more appropriate to describe the 
implementations discussed in the next section. 

In the Abelian case, the system with nontrivial topol- 
ogy, i.e., with holes, acquires a degenerate ground state, 
since a fluxon located inside the holes has a vanish- 
ing energy. Such a fluxon is created by the operator 
ri(Qb)(E7 1~ab\^) where the contour 7 in A connects the 
inner and the outer boundaries. This operator is also 
equal to ri(ab)e7 ex P(^ 7r ^ afc / rt ) an< ^ ^ s conserve d by the 
Hamiltonian (|14|l . Its eigenstates correspond to linear 
combinations of various degenerate ground states with 
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fluxons present or absent in the holes. This applies di- 
rectly to the non-Abelian model of the present section. 
The degenerate ground states are therefore identified by 
the value of II(ab)67 ^ ab ' ^ n s ^ e re P resen t a tion, it is 
just the group element ging^ut where jm, g ou t represent 
the states of the inner and outer boundaries respectively. 
These arguments hold provided that each boundary is 
described by a single group element g i.e., that all sites 
belonging to it are connected by nonfluctuating bonds. 
In the Abelian case discussed in Sec. IIIII the boundary 
(denoted by the dashed line in Fig. 0) is just a supercon- 
ducting wire. 

The same arguments as the ones developed in the 
Abelian case |7j show that only nonlocal operators in g a b 
result in the decoherence of the topologically degenerate 
states. This shows that these degeneracies are protected 
from decoherence provided that no physical operator cou- 
ples to the state g a of an individual site. 

V. IMPLEMENTATIONS 

As explained in Sec. IIIII the array that implements 
directly (in a natural basis) the gauge-invariance con- 
straint for the Abelian 2„ group is physically an insula- 
tor while the superconducting array becomes equivalent 
to the same gauge theory only after a dual transforma- 
tion. Alternatively, the superconducting array can be 
thought of as implementing directly the modified Gauss 
law discussed in Sec IV B. In the non-Abelian case these 
two theories are no more dual to each other, so we discuss 
these cases separately. 

The common element of all arrays discussed here is 
the Josephson-j unction circuit that is characterized by 
the phase difference A(j> which can take (in a classical 
limit) n > 2 discrete values A<p = 2-Km/n] quantum pro- 
cesses leading to transitions between these states. The 
actual construction of this element for n = 3 (that we 
shall call masu |23J below) is shown in Fig. ^| It is built 
from n more basic elements consisting of two Joseph- 
son junctions connected in series. The Josephson en- 
ergy of this element (as a function of the phase differ- 
ence 5cj> between the ends) contains all higher harmonics: 
Vo(5(f>) = —2Ej\cos8<fr/2\. Now, connecting these ba- 
sic elements in parallel as shown in Fig. 1 101 and applying 
the magnetic field that shifts the relative phases of neigh- 
boring circuits by 2ir jn gives the desired 27r jn periodic 
potential. The nonzero charging energy of the junctions 
leads to transitions between these classical states whose 
amplitudes are calculated in Appendix A. 

A. Gauge-invariance implementation 

Here we limit ourselves to the gauge theory with the 
simplest non-Abelian group D% which is also the permu- 
tation group of three elements S3. The array implement- 
ing this model is shown in Fig. \7\ We begin by explaining 



the qualitative reasoning that leads to this complicated 
construction. In the S3 gauge theory the state of each 
link has to be described by the permutation of three el- 
ements. To implement this in a Josephson-junction ar- 
ray, we construct each bond from two elements (masus) 
described by discrete phases differences A(j> a = 2irn/3 
(a = 1,2). Such a bond can be viewed as a discrete 
transformation from the state of one site characterized 
by two discrete 2irn/3 phases to the state of another site; 
it is described by one of the nine elements of the Z3 x Z3 
group. In order to eliminate three extra states and to 
get the permutation group, we need to introduce a con- 
straint requiring that two phases of each site be different. 
This is achieved if each site of the lattice contains an ele- 
ment that excludes the configurations with equal phases. 
The simplest of such an element is a triangle made of 
three Josephson junctions frustrated by a magnetic flux 
$o/2. The ground state of this triangle corresponds to 
the phases <j) a = 2im a /3 with different n a . Thus, the 
state of two masus attached on both sides to the frus- 
trated triangles [as shown in Fig. Efb)] is in a one to one 
correspondence with the elements of S3 group. 

The next step in the construction of the array is 
to ensure that the low-energy states satisfy the gauge- 
invariance constraint. This is true if the Hamiltonian 
is dominated by a large "kinetic" term which exchanges 
two discrete phases on each triangle because such a pro- 
cess changes the states of all three bonds connected to 
the triangle, further, it is equivalent to the action of a 
transposition on the states of these three bonds 

2 

H T = ~r Y,Y,I[ T ^ k )> ( 48 ) 

a fc=0 6(a) 

where ro is the tunneling amplitude of this process: 

ro^ErE^^i-So^Y (49) 

where So = 6.10 if junctions in triangles have the same 
Josephson and charging energies as the ones in masus (see 
Appendix B) . Note that this numerical value is sensitive 
to the number of nearest neighbors of the effective lattice 
and is minimized for the hexagonal lattice. As will be 
clear from the following, in order to be described by the 
gauge theory, r should be sufficiently large which is a 
case easier to achieve for the hexagonal lattice. 

The wave functions minimizing the kinetic energy 148|) 
are the ones that satisfy the gauge-invariance constraint 

Y[T ab (rt k m = \y), (50) 

b(a) 

since local transpositions generate the whole S3 gauge 
group. The array is faithfully described by the gauge 
theory if ro sets the largest energy scale in this problem. 
This is not satisfied if the bonds of the array are just 
the ones shown in Fig. E[b), because in this case the 
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FIG. 7: (a) The insulating array implementing non-Abelian 
insulator. Each bond of the hexagonal array is built from 
two circuits, each consisting of one masu (rhombus with di- 
agonal) and one weak link (shown as a Josephson-j unction 
loop with two junctions). The sites of the array are formed 
by triangles frustrated by a magnetic flux $o/2. (b) A sim- 
pler construction of one bond which states are in one to one 
correspondence with the elements of the permutation group. 
A and B denote the leads that can be used to measure the 
array properties in the experiments similar to those discussed 
in Ref. 



loop is ej = Ej($/<&o — 1/2) and its charging energy is 
€c = Ec/2- The effect of these small energies on the tun- 
neling process that corresponds to transposition is weak 
and can be ignored. The presence of weak links decreases 
the energy gap between the six states of the bond asso- 
ciated with the permutation group and the remaining 
three states that were eliminated from the low-energy 
spectrum. However, as we show below, we can choose 
the energies of the weak links, ej, ec, so that this gap re- 
mains large, but the amplitude of the tunneling process, 
ro, remains larger than other scales in the low-energy the- 
ory so that the leading term in the effective Hamiltonian 
is indeed given by Eq. 

The Hilbert space satisfying this constraint is still 
huge: in the absence of other contributions to the Hamil- 
tonian (generated by weak links) all states satisfying the 
gauge-invariance constraint are degenerate. This degen- 
eracy is lifted when these subdominant terms are taken 
into account. Weak links give two types of terms: they 
lead to the dynamics of the discrete variables describing 
the state of each link and to diagonal potential-energy 
terms. Let us start with the latter. The main contri- 
bution comes in the second order in ej/ec, it is given 
by 



U({ui}) 



2e c 



2tt 



(51) 



where Ui are discrete variables describing the phase dif- 
ferences on elementary masus and along the sides of the 
triangles. This energy is minimized if this sum is mo- 
dulo 3, i.e., if the state of the links is described by the 
permutation group "acting" on the state of triangles. We 
assume that the scale of this interaction is relatively large 

and comparable to ro, i.e., that > ro. This interac- 
tion does not lift the degeneracy of the Hilbert space of 
the discrete variables that are described by elements of 
the permutation group. This degeneracy is only lifted 
by weaker processes which involve longer loops of weak 
links around each hexagon. Such a loop contains six weak 
links, it results in the contribution to the potential energy 
similar to the first term in Eq. (|26|l : 



product of permutations around the elementary hexagon 
is unity and it is the energy that is needed to violate this 
constraint that sets the largest energy scale in the array. 

In order to decrease the energy scale associated with 
the Ilhez 9 ab 1 configurations, we need to introduce 
weak links (shown as loops containing two Josephson 
junctions in Fig. |7|) characterized by small Josephson and 
charging energies ej <C ec <C Ec ^ Ej. In order to avoid 
a mixture of Josephson junctions with different param- 
eters (Ej,Ec) in the same array, we assume that these 
weak links are constructed from the two usual Joseph- 
son junctions with a large Josephson energy Ej > Ej 
and correspondingly small charging energy Ec ^$ Ec 
connected in parallel to form a loop frustrated by the 
magnetic flux. The effective Josephson energy of such a 



V = ^ V (9ab9bcg C dgdegefgfa)- (52) 
abcdef 

To calculate the value of the interaction V(g), we note 
that physically this energy is due to the energy of the 
Cooper pair that starts from one corner of the triangle 
and moves around the hexagon across the weak links. 
Furthermore, if g is a transposition, the sum of the dis- 
crete phases along contours beginning at two corners of 
the triangle is equal to ±27r/3 and one such a sum is 0, 
while if g is a cyclic permutation, all these sums are equal 
to ±2tt/3, so V(e) = 0, V{r( k ) = 2v /3, V{( k ) = v , 
where coefficient u is the strength of the effective Joseph- 
son contact formed by two parallel loops of six weak junc- 
tions with a phase difference 27r/3 across. To calculate 
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the value of vq we note that two parallel loops are equiv- 
alent to one loop with effective Josephson and charging 



eff 

energies ej J 



2c 



J, e 



eff 



ec/2 and use the result of 



the Appendix C for the energy of one loop 



20 I e c 



(53) 



We now turn to the dynamical processes which involve 
phase slips of weak links. These processes change the 
state of each link individually, similar to the second term 
in Eq. (|26l) . The potential energy involved in such a 
process is dominated by the energy of each masu while 
the charging energy is due to the weak junction (because 
its charging energy is smaller, it dominates the dynamic 
term in the action). Thus, the amplitude of such a pro- 
cess is 



For small ec this amplitude becomes exponentially small 
and so these terms can be neglected in the effective 
Hamiltonian J2fjJ. 

Thus, we conclude that the array shown in Fig. [7| 
provides a realization of the non-Abelian gauge theory 
discussed in Sec. |W] if ej -C e c "C E c ^ Ej and 



2e c 



r > 



20 



3ejV 
iec) 



(55) 



B. Kinematical Gauss law implementation 

As we argue below, the straightforward implementa- 
tion of the kinematical non-Abelian Gauss law is not 
possible because it does not protect the system from the 
coupling to the environment. In this respect, the non- 
Abelian theories are crucially different from the Abelian 
ones. In order to explain this difference, we first recollect 
the physical arguments that show that Abelian arrays 
(similar to the ones suggested in Ref. Q) are optimally 
protected from the environment and then show why non- 
Abelian ones are different. 

For Abelian theories, the Gauss law (0J and (|35|l is even 
easier to implement in the physical system than a gauge 
theory. As discussed in Sec. II, in order to implement 
an Abelian group Z„ on the triangular lattice shown in 
Fig. 1 one might assign elements (masus discussed in the 
preceding subsection) that are 1i:jn periodic in super- 
conducting phase to the links of the dual hexagonal lat- 
tice. In this system the 'site' representation acquires a 
transparent physical meaning - the state of each island 
is described by the superconducting phase of each site of 
the hexagonal lattice. Different topological sectors of the 
system with hole(s) correspond to different phases of the 
inner boundary. The existence of the protected subspace 
is related to the U(l) physical invariance: superconduct- 
ing phase of each island (and of the inner boundary) has 



no meaning by itself, so only the difference of the phases 
of the neighboring sites might be coupled to physical op- 
erators. 

The operator conjugated to the superconducting phase 
is the physical charge of the whole island and this 
presents some danger since it can be coupled to the exter- 
nal electric fields, especially if it results in a direct cou- 
pling of the state of the inner boundary to the electrical 
fields. Such a coupling, however, becomes exponentially 
small in the system with a long range order in super- 
conducting phase, more precisely in the systems where 
(exp(zn0)) is close to 1. In these systems the total charge 
of the boundary fluctuates strongly, while the conjugate 
operators appearing in the discrete gauge theory corre- 
spond to the total charge modulo n and thus are decou- 
pled from the physical electric fields. We compute their 
weak residual effect on the degeneracy of the protected 
states in Appendix D and find that it quickly becomes 
exponentially small [see Eq. i|D6l) ] with the number of 
links Z connecting the boundary to the inside of the ar- 
ray, i.e., in the thermodynamic limit. 

As mentioned in the preceding section, the most 
straightforward way to generalize this Z n construction 
to the simplest non-Abelian group S3 would be to con- 
struct the bonds characterized by the elements of the 
permutation group and the sites that can be viewed as 
results of these permutations. This is achieved in the ar- 
ray shown in Fig. [S] where bonds are constructed from 
pairs of masus and the sites from triangles frustrated by 
the magnetic field. This array differs from the one shown 
in Fig. 0by the absence of the weak links. All classical 
states of this array are uniquely specified by the configu- 
rations of permutation group elements which satisfy the 
constraint Yihex 9 ab = e. A small but nonzero charg- 
ing energy of each contact would lead to the quantum 
transitions between different classical states as described 
by the Hamiltonian l|3()|) . In this physical system, these 
transitions correspond to the tunneling of the phases on 
each triangle from one configuration to another. 

Thus, in ideal conditions, this system is described by 
the kinematic Gauss law model of preceding section. The 
situation is not so good for a realistic system because in 
such implementation, one can not exclude the appearance 
of physical operators coupled directly to the state of each 
triangle. In particular, this could distinguish between 
states with different chiralitics (i.e., in which the phase 
winds clockwise or counterclockwise around the triangle). 
This is bad news for the effective gauge field model since 
the phase differences on distant triangles are related by 



A0i = Acj> 2 + ysign fjjsab^ 



(56) 



where sign denotes the sign of the permutation and in- 
dices 1 and 2 denote two distant triangles. Thus, the 
terms in the energy that depend on the phase difference 
Acb results in nonlocal terms in the effective Hamiltonian 
when expressed as a function of the gauge fields g a b- Such 
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the noise 





(b) 



(a) 



FIG. 8: (a) Implementation of the non-Abelian theory with 
group D3. The elementary block contains two elements 
(called masus and shown as rhombi), each element is peri- 
odic in the phase across it with periodicity 2ir/3, and one 
triangle which contains Josephson junctions in each side and 
is subjected to magnetic flux $o/2. This ensures that two 
phases on the ends of each block are different. The construc- 
tion of each masu is shown separately in (b) . The periodicity 
is achieved applying magnetic field with flux $0/3 through 
each elementary triangle. 



terms ruin the topological protection, especially if these 
triangles belong to inner and outer borders. To be more 
specific we shall discuss below the effect of 



V dl . 



F d cos(A0 1 - A0 2 ), 



(57) 



which is actually the most dangerous term of this type 
since it changes the relative energies of different topolog- 
ical sectors of the array shown in Fig. |H1 by an amount 
proportional to Vd- 

The effect of V d is. on the array with weak links (dis- 
cussed in Sec. IV A|) is exponentially weak because in this 
system the relation between phases also contains a con- 
tribution from the continuous phases across weak junc- 
tions which, in the absence of the superconducting long 
range order, fluctuate strongly and eliminate the corre- 
lation between A^i — A</>2 and \\c 9ab- Indeed, in this 
case, the relation between A<^i and A</>2 acquires an ad- 
ditional term which is due the continuous part of the 
phase across the weak links. 

- A0 2 = ^sign(Q g ab ) + ^ 50, (58) 



where 5<f> is the difference of continuous phases across two 
weak links constituting one bond of the effective lattice. 
In the insulating limit of this model (ec S> £j) the contin- 
uous phases fluctuate practically independently on each 
bond with (cos(6<f)) ~ (cj/ec), suppressing the effect of 



V. 



eff 



dis. 



V d 



2tt 



sign Yl 9ab 



(59) 



where L is the length of the contour connecting triangles 
1 and 2. 

Alternatively, one can argue that the mapping of the 
physical system onto the effective model with variables 
defined on the bonds of the lattice does not produce non- 
local (dangerous) operators if the states of the site are 
indistinguishable. In the insulating array, this is achieved 
by making the phase of each triangle (effective site) fluc- 
tuate strongly. In the Abelian Z n array, this is true 
because absolute values of the phase (which distinguish 
between different states of a given site) are not measur- 
able, so generally a system in which each site is charac- 
terized only by a single phase is a good candidate for an 
array that can be mapped to a bond model. The problem 
is that the straightforward implementations of this idea 
lead to a 2 n group. 

In order to add a reflection to the Z„ group, one needs 
states that are characterized by both the discrete value 
of the phase and the sign of the conjugate variable, q. 
This is possible to achieve if one starts for the system 
where n x m discrete phases are allowed, and include a 
large interaction that forms a degenerate doublet of chiral 
states from the subset of m discrete states. To be more 
specific, consider the low-energy states that are described 
by the wave functions 



# ± (Z) = exp(±2Trik/m)\l,k), 



(60) 



k=0 



for I < n. Because the reflection changes the chirality of 
each state, it is not reducible to a simple rotation, so that 
combined with the phase rotations, it gives the funda- 
mental representation of D n group. We now discuss the 
construction of the simplest array that has these general 
features. The minimal value of m that allows one to dis- 
tinguish left and right states is m = 3. The minimal value 
of n that allows a non-Abelian group is n — 3, however, 
n = m is very difficult to implement because it can not 
be achieved connecting in series the elements with n and 
m periodicity (that would give n = m periodicity instead 
of n 2 one). Thus, the simplest array of this type corre- 
sponds to n = 4 and group -D4. A possible realization is 
shown in Fig. Here each link consists of two elements 
which have rather different values of Ej/Ec- the masu 
that is 27r/4 periodic operates in the phase regime with 
Eji/EcA > 1 while the masu that is 27r/3 periodic has 
a different relation Ej^/Ecd, ~ 1- More precisely, the 
amplitude mixing states that differ by 27r/4 is 



^ 4 ^ / 4 4 exp(-3S 4 ), 



(61) 



where S4 = yj Ej4 / Eca must be much smaller than the 
energy gap 



A«£%V 3 4 exp(-3S 3 ) 



(62) 
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between the doublet of the chiral states and the re main- 
ing p — state of the m = 3 set (S3 = yEjz/Ecl). 




(a) (c) 

FIG. 9: (a) Robust implementation of the non-Abelian 
group D4 involving elements in phase and charge regime, (b) 
Rhombi denote elements that are 27r/3 periodic and whose 
charging energy is relatively high, (c) Big squares denote ele- 
ments that are 2-7r/4 periodic in phase, the junctions of these 
elements are in a quasiclassical regime. 

It is crucial for this design that the states with p = ±1 
[here p has the same meaning as in Sec. Ill, e.g., in Eq. 
(|ll)fl ] have equal energies which are less than the energy 
of the state with p = 0. This happens if the transition 
amplitude between the classical minima is real and neg- 
ative. Generally, the amplitude of the transition that 
increments the superconducting phase by 27r/3 contains 
the phase factor 2ttQ /3 where Q is the average charge of 
the superconductor whose phase has been changed (here 
and below we measure the charge in the units of Cooper 
pairs). We require that this phase is (7r + 27riV) where N 
is arbitrary integer. This implies that the system should 
be biased so that the average charge of the island that 
experiences this transition is Q = 3/2 + 3N. 

We now discuss the effect of the induced charge on the 
amplitude of the transition between I = 0, . . . , 3 discrete 
states. For simplicity, we shall assume that the islands 
which experience these transitions coincide with the is- 
lands whose phases change by 27r/3 and, thus, have the 
same charge. These transitions increment phase by 7r/2; 
so that their amplitude acquires a phase (3/2 + 3iV)7r/2; 
in the notations of Sec. IIVBI = exp(iQ7r/2)|r|. This 
amplitude has positive real part for Q = 9/2, 15/2, . . .. 
The transition between states with p = ±1 appears nat- 
urally if the energies of the classical minima with differ- 
ent k = 0, 1, 2 are not exactly equal because, in the mo- 
mentum representation, this induces a transition between 
these states. Such a difference is induced by the devia- 
tions, (5$3 = $3 — <E>o/3 of the actual flux through the 
elementary rhombus from its ideal value, $0/3. Thus, 
in these conditions the low-energy dynamics of such a 
system is described by exactly (|37f) with the parameters 
TV, and r T = Ej^S^^/^o. We emphasize again that the 
low-energy sector has to be separated from higher-energy 



excitations by a large gap A ^> r^, r T . 



VI. CONCLUSION 

In this work, we have proposed several possible im- 
plementations of lattice models based on discrete non- 
Abelian local symmetry groups. The first key ingredient 
in all these constructions is the existence of a large but 
discrete manifold of equivalent "classical" states. This 
manifold appears if the system is built from heavily 
frustrated elementary blocks, such as Josephson-junction 
loops threaded by a well adjusted magnetic flux. Fur- 
ther, one needs to select an interaction between these 
blocks that imposes a set of local constraints in the low- 
energy sector. We have found two inequivalent (but 
somewhat similar) types of constraints compatible with 
the non-Abelian symmetry group. The first type of con- 
straint corresponds to the dynamical selection of the 
states which are invariant under the gauge group trans- 
formations. This is most naturally implemented in in- 
sulating arrays where it can be made equivalent to the 
selection of states that minimize the kinetic energy as- 
sociated with the local tunneling processes between clas- 
sically equivalent states. The second type of constraint 
emerges more naturally in superconducting Josephson- 
junction arrays (or the ones with small phase fluctua- 
tions), where group elements associated with the bonds 
can be viewed as transformations between the states of 
the sites of the array. Because the state of each site is 
uniquely defined, the product of group elements around 
a closed loop is the unit element. We have called this 
condition a kinematical Gauss law. This constraint be- 
comes equivalent to the local gauge-invariance condition 
in Abelian theories where a duality transformation maps 
one onto the other (Sec. IIII|) but for a non-Abelian sym- 
metry it leads to a different theory in which the roles of 
charges and fluxons have been reversed (Sec. II V [I . 

In both types of models that we considered, the ele- 
mentary excitations within the constrained subspace are 
local in space and have a finite-energy gap. We call these 
excitations "fluxons". Generally, these fluxons can be 
trapped inside a larger hole prepared in the array lead- 
ing to a degenerate ground state in the thermodynamic 
limit. In a finite system, this degeneracy is lifted by an 
exponentially weak tunneling escape process. It is also 
robust with respect to all noises represented by local op- 
erators such as local electric fields (generated by stray 
charges) or stray magnetic fluxes: the matrix elements 
of all these operators are exponentially suppressed. This 
remarkable property may open the way to the interesting 
applications such as storage of quantum information for 
which a low decoherence rate is a must. We emphasize 
that the physical nature of these fluxons depends, how- 
ever, on the details of implementation and the type of 
the constraint: the canonical gauge invariance is more 
easily implemented in insulating arrays in which fluxons 
have the physical meaning of fractional vortices while 
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kinematical Gauss law is more naturally implemented in 
superconducting arrays in which fluxons carry physical 
charges. 

Besides fluxons, in any physical array, it is possible to 
build excitations that violate a local constraint (of course, 
at some energy cost); in the gauge theory language these 
excitations are pseudocharges. In the canonical gauge 
models pseudocharges may provide a multidimensional 
representation of the discrete gauge group while fluxons 
are labeled by their conjugacy classes. In the kinematical 
Gauss law models the role of fluxons and pseudocharges 
is reversed. The existence of these higher-dimensional 
representations leads to a very interesting statistical in- 
teraction between fluxons and pseudocharges which may 
be utilized for making exact transformations needed in 
quantum computations. 

Many interesting questions remain to be addressed, we 
shall list just a few here. First, it is not clear whether it is 
absolutely unavoidable to emulate gauge theory in order 
to get the topological protection and non-Abelian statis- 
tics. What is the most general class of theories that has 
these properties? Second, is it possible that some sim- 
pler overfrustrated arrays "spontaneously" develop the 
gauge invariance in the low-energy sector? |24j . Such an 
array would be more experimentally accessible and eas- 
ier to control. Coming back to the proposed arrays it 
would be very interesting to investigate the experimen- 
tal manifestations of the nonlocal nature of the fluxon 
quantum numbers which become apparent only when a 
few of them are present. As discussed extensively in the 
mathematical literature |!8lll9| . these quantum numbers 
interpolate between a set of conjugacy classes (one for 
each fluxon) and a set of group elements (which is not 
a gauge- invariant object by itself). It would be impor- 
tant to realize how these subtle but beautiful concepts 
can be revealed in a physical experiment, and to design 
some setup that could verify the non-Abelian statistics 
of the excitations in these models. Finally, it would be 
very interesting to adapt the construction of the Chern- 
Simons term (which binds fluxons and charges) known for 
a continuous groups to the discrete finite gauge groups 
(even Abelian ones) discussed here and to provide a cor- 
responding physical implementation. 
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bus with diagonal that is placed in the magnetic field with 
the flux 2$o/3 through the whole rhombus. In the follow- 
ing we shall refer to this object as masu. As a function 
of the phase difference ip across the diagonal the total 
Josephson energy of such masu is: 

E(ip) - -IE j [|cos(^/2)| + |cos(V>/2 + 7r/3)| 

+ |cos(^/2-tt/3)|]. (Al) 

It is periodic in ip with the period 27r/3 and has three 
classically degenerate minima: ip = 0, ±27r/3. 

The amplitude of transitions r m between these minima 
can be computed in a quasiclassical approximation, as 
r m w Ej Eq exp(— So), and is equal for all transitions. 
To calculate the action So, we choose two independent 
superconducting phases as shown in Fig. 1101 and rescale 
t into t/y/EjEc- In terms of these variables the action 
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2a/3cos(6>) 



dt. 



(A2) 



The minimum of this action corresponds to the quasi- 
classical equations of motion: 



d^d 8 . 8 . „. . ,. , , ... 

dfi = gcos^cos^), (A3) 
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8sin(0) sinf». 



(A4) 



Solving these equations numerically, we get the expres- 
sion for the transition amplitude in one masu 



Ej /4 ^ /4 exp(-S'o) So = 1.22 



(A5) 
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APPENDIX A: TUNNELING IN A SINGLE 
RHOMBUS 

Here, we calculate the transition amplitude between 
the degenerate classical states appearing in a single rhom- 



FIG. 10: Tunneling in the elementary building block of non- 
Abelian arrays. In the quasiclassical approximation, tunnel- 
ing is described by the simultaneous evolution of all phases 
that are parametrized by two independent variables, 8 and (p. 
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APPENDIX B: 



TUNNELING OF FOUR RHOMBI 
AND TRIANGLE 



Here, we calculate the transition amplitude of four ma- 
sus connected by a triangle as shown in Fig. ^2 If the 
junctions constituting the triangle are weak, ej <C Ej 
this amplitude can be obtained in the perturbation the- 
ory as a product of two transitions that happen consec- 
utively: one involving the left pair of masus and another 
the right one. The amplitude for each transition is re- 
lated to the amplitude of the transition in a single masu 
(|A5|I because the quasiclassical equations of motion for 
the variables describing two masus joined at the vertex 
are the same as for a single masu. Thus, this tunneling 
amplitude contains twice the action (|A5(I : 

r' m « Ey A E l J" exp(-5 ) 5 = 2M^. (Bl) 

In this limit, the amplitude of the full process is 

(r'm) 2 



ro 



ej 



(B2) 



where again ej denote the weak Josephson couplings in- 
side the triangles. In a general case of triangle junctions 
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-2cos(26>) + cos(46>)] 



(B3) 



In the most natural case the junctions of the triangle 
have the same Josephson current per area and the same 
capacitance per area as the junctions of the masu, so 
ej/Ej = Ecj&c = e. Minimizing the action, we get the 
equations of motion 



2 2 
—= sin(0) — - cos(0) cos(0) 
V 3 3 



[sin(40) - sin(20)] 



-jZ = 8sm(0)sin(». 



(B4) 
(B5) 



Solving these equations numerically, we obtain the di- 
mensionless action 5(e) as a function of e shown in Fig. 
II 21 which should be used in 



£ 3 /X /4 exp 




(B6) 



to obtain the value of the transition amplitude. We 
observe that the action 5(e) varies only by 25% [from 
5(0) = 4.91 to 5(1) = 6.10] when e is increased from 
to 1, so the amplitude of the transition is not too sensitive 
to the strength of Josephson junctions in the triangles. 



FIG. 11: Basic tunneling process in the full lattice. 




0.2 



comparable with the junctions constituting masus, we 
need to add these terms to the action and solve the full 
set of equations of motion in order to obtain the action 
corresponding to this process. The solution of these equa- 
tions describes a simultaneous variation of phases in all 
four masus, so the most general action has the form 



5 



1 



2^3008(0) 



-2cos(</>) sin(0) 



ZE C 
8e c 



FIG. 12: Dimensionless action of the tunneling process. 



APPENDIX C: EFFECTIVE JOSEPHSON 
COUPLING OF A CHAIN 

Here we calculate the effective coupling by a chain of 
(ri + 1) weak Josephson contacts. Denoting the phase 
across each contact by tp^, the phase across the last con- 
tact is 4> — ifi. After a time rescaling the system is 



19 



described by a dimensionless Lagrangian 
L = Lq + Lj, 



Lj = A 



cos Lp l + cos I <t> - ip % 



(Cl) 
(C2) 

(C3) 



where A = ej/(8ec)- Let us assume that A <C 1 and 
treat Josephson energy in perturbation theory. The bare 
Hamiltonian is quadratic in charges, conjugated to the 
phases ipi 



l \ 2 



(C4) 



The ground state of this Hamiltonian corresponds to all 
charges = while in the low-energy states some qi = 
±1. The energy of the state with k charges (/, = 1 is 



e fc = fc(n + l-fc)/[2(n+l)]. 



(C5) 



To find the effective coupling provided by this chain we 
need to calculate 



in the lowest order of the perturbation theory in 



Vim = 3 E eX P(*V*) + h ' C - 

i 

In the charge representation, one has 
W(4>) = -\(l[q+ 



COS( 



(C6) 



(C7) 



(C8) 



where qf is the operator increasing the charge of junction 
i by 1 and 



V ir 



E< 



/i.e. 



(C9) 



In the leading order, we can limit ourselves to the com- 
ponents of the ground-state wave functions that contain 
charges —1 at some k junctions and (n — k) opposite 
charges at the remaining junctions. These components 
have coefficients that can be found from a recursive 
relation 



so we get 



k\ 

*fe = t^— *fc-i. 

2£fc 



(CIO) 



(Cll) 



Expressing the average (J\ i q~l) through these compo- 
nents via 



ru + =e 



fc=0 



{n-k)\k\ 



k ^ n — ki 



(C12) 



we get 



(C13) 



Converting back to unrescaled variables we conclude that 
the effective coupling between the ends of the line con- 
sisting of n junctions is 



W n (</>) 



n / \ n— 1 



APPENDIX D: EFFECT OF THE ELECTRIC 
NOISE ON THE PROTECTED SUBSPACE IN Z„ 
ARRAYS 

Quantitatively, consider the island coupled to Z links 
by 2tt /n periodic elements. Each element is characterized 
by the potential energy V e i{4>) — J2k Vo(4>+2i:k / n) which 
has minima at <f> m — 27t(to + \)/n. Near each minima 
the potential is approximately Gaussian: 



V el ^-Ej{cj>-<P m f, 



where 



1 71—1 



k=0 









cos 













(Dl) 



(D2) 



is a numerical coefficient that is equal to \j\pl for n = 
2 and n/ir for large n. The dynamics of small phase 
fluctuations is described by the Lagrangian 



E 

m 



32E C 



{<t>i-4>jf + V e l(<j>i 



(D3) 



where Ec = e 2 / (2C) is the charging energy of each junc- 
tion. In this approximation, each phase difference is an 
independent Gaussian variable which is described by the 
wave functions 



= cxp(- T 



ncEj 



2 )- 



(D4) 



The canonically conjugate charge fluc tuations arc the rc- 
fore also Gaussian with width J (ncEj)/(l6Ec)/2; 
charge fluctuations on different junctions add up to the 
total charge fluctuation on the island with the Gaussian 
width 



((Q-QY) = t 



Z / ncEj 



8 V Ec 



(D5) 
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If these fluctuations are restricted to a given value mod- 
ulo n, they give slightly different results. Consider, for 
instance, the difference between average square of the 
charge between the sectors with different Q modulo n 

W 2 ) fc -W 2 )^exp^y^). (D6) 

This difference quickly becomes exponentially small for a 
large number of links connecting a given island. Because 



a topologically protected space can be viewed as spanned 
by different states of the inner boundary of the array 
connected to the interior by Z 1 links, the electrical 
fields have exponentially weak effects on the topologically 
protected subspace. 
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